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Abstract. Tropical Kelvin waves have been observed previously in ozone mixing ratio 
data from the SBUV (Solar Backscatter Ultraviolet) and LIMS (Limb Infrared Monitor 
of the Stratosphere) instruments on board the Nimbus-7 satellite. The present study 
investigates Kelvin wave features in total column ozone, using version 6 data from the Total 
Ozone Mapping Spectrometer (TOMS) instrument (also on board Nimbus-7). Results show 
eastward-propagating zonal waves 1-2 with periods ~5-15 days, amplitudes 3-5 DU, and 
latitudinal symmetry typical of Kelvin waves. A simplified model calculation suggests 
that the primary source for the perturbations is slow Kelvin waves in the lower-to-middle 
stratosphere. Maximum Kelvin wave signatures occur in conjunction with westward lower- 
to-middle stratospheric equatorial zonal winds [a quasi-biennial oscillation (QBO) wind 
modulation effect]. 
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Introduction 

Many observational studies have been done on tropical Kelvin waves, with some of the 
earliest being Wallace and Kousky (1968), Maruyama (1969), and Angell et al.(1973). The 
present investigation examines TOMS total column ozone for signatures of these waves. 
As will be seen, the lower stratosphere appears to be the primary region responsible for 
Kelvin waves in total ozone. 

Stratospheric Kelvin waves in ozone mixing ratio were studied by Prata (1990) using 
SBUV data for January 1979. That study reported upper stratospheric Kelvin waves in 
zonal waves 1 and 2 with periods ~1 week and vertical wavelengths ~ 15 km. A more 
recent study of Kelvin ozone waves by Randel and Gille (1991) examined LIMS data for 
January-February 1979 and SBUV data for January 1979— February 1986. Typical periods 
of upper stratospheric Kelvin waves were found to vary from 5 to 15 days with vertical 
wavelengths comparable to those found by Prata. Randel (1990), using LIMS data, showed 
that there are two primary regions of Kelvin waves in ozone mixing ratio: (1) the upper 
stratosphere (where photochemistry dominates), and (2) the lower stratosphere (where 
transport dominates). In the middle stratosphere (~10 hPa), Kelvin wave features axe 
virtually nonexist because of vanishing vertical gradients in ozone mixing ratio. 
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For investigating interannual variability of Kelvin wave signals in the lower strato- 
sphere (where vertical wavelengths ~6-10 km are common), there axe problems with using 
the SBUV data. As noted by Randel and Gille, vertical resolution is ■ — 8-10 km, and there 
is also an apparent effect by the processing algorithms that restricts detected vertical wave- 
lengths to 27 km or larger. In contrast, the TOMS data, as will be seen, appear to be 
sensitive to Kelvin waves in the lower stratosphere. 

Three 60-day case studies will be presented below, with each exhibiting features typical 
of tropical Kelvin waves. Two episodes are near solstice, beginning in July, and the third 
near equinox, beginning in March. 


Data and analysis 

This study uses TOMS version 6 daily data from the NASA/Goddard Space Flight 
Center made available on CD-ROM disk in a format of l°latitude by 1.25°longitude bins, 
covering 1 November 1978 to 31 January 1992. These data were regridded in this study 
onto a 5° x 5°mesh extending from 85°N to 85°S in latitude and 180°to 175° W in longitude. 
Further details concerning regridded TOMS data and spectral methods are discussed by 
Stanford and Ziemke (1993). 

Monthly mean zonal winds from Singapore (1.37°N, 103.92°E) [B. Naujokat, Free 
University of Berlin, personal communication] for years 1979-1991 at pressure levels 70, 
50, 40, 30, 20, 15, and 10 hPa were used in this study. 


Kelvin wave modeled amplitudes in total ozone 

A model for tropical Kelvin waves in total column ozone is formulated using linearized 
forms of the thermodynamic and tracer continuity equations appropriate for Kelvin waves. 
The two equations are, respectively, 


T[ + xTT' x + w'HN 2 /R = -TV" 1 (2) 

x't + ux'« + w'x, = -iY - er (3) 

where T is absolute temperature, \ is volume mixing ratio, H=scale height (7 km here), 
R=gas constant (287 J-kg -1 K -1 ), JV= buoyancy frequency (taken to be 2 x 10 -2 s -1 ), r -1 
is a radiative damping coefficient, T and 6 are linear photochemical production/destruction 
parameters, the over bar denotes a zonal average, and subscripts denote partial differenti- 
ation. All other symbols are standard. Zoned mean wind u is a constant, but remains as 
a variable parameter in this model. 

In (2) we use an approximate empirical form for the radiative damping coefficient given 
by r- 1 (units: s" 1 ) = A( 1 - e”®/*.).^/^) 1 - 8 , where A = 3.7 x 10" 6 s~\ £=14.8 (km), 
z is altitude (km), and X z =vertical wavelength (km). This expression was constructed 
upon examining the diagrams of r -1 in Fels (1982, 1984). 

In (3), empirical forms for 9 and T are also used, given by 6 = 3.33 x 10 -16 • ( z — 40), 
and T = 1.67 x 10“ 8 • (z — 40), respectively. Here, z is altitude in units km, and units for 9 
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(r) are K -1 s -1 (s -1 ). Both 6 and T are taken to be zero for z < 40 km. These expressions 
are simple linear approximations to the 9 and T shown in Fig. 3 of Hartmann and Garcia 
(1979). Standard interpolated UMKEHR profiles (R. McPeters, NASA/GSFC, personal 
communication) for 250 DU in the tropics were used to estimate x z i n (3)- The bottom of 
UMKEHR level j (j — 0, 1, 2, 10) is given by pressure (1013.25)-2 - -' hPa. 

Perturbations in total column ozone (units: meters) Eire given by 

/ OO y-OO 

dz • p ■ x' — * / dz-e-*' H -x' (4), 

o Jo 

where p , the total mass density, was approximated in (4) by 1.3 • e~ z / H (units: kg-m -3 ). 

For induced Kelvin wave perturbations, the form of x' is taken to be the real part of 
x'{ x iyi z ,t) = x(y,z) • e z / 2H ■ where the exponential terms were chosen to 

balance those of linear Kelvin wave temperature perturbation solutions (see, for example, 
Andrews et al., 1987). Based on observations, Kelvin wave temperature disturbances in 
this model were all taken to be 2.8 C°at 4 scale heights. 

From (4), peak perturbation amplitudes in total column ozone are found to be ~l-2 
DU for eastward propagating linear Kelvin waves 1 or 2, with vertical wavelengths 8-32 
km (see Fig.l). This linear model appears to break down for short vertical wavelengths; 
only wavelengths longer that 8 km are plotted in the figure. According to Bhartia et al. 
(1984), local errors in unfiltered TOMS data are ~l-2%. In comparison with the model, 
upper stratospheric Kelvin waves with long vertical wavelengths (> 15 km) would induce 
amplitudes too small to be detected in TOMS data. 


Kelvin wave signatures in TOMS 

TOMS power spectra for fast (periods < 15 days) zonal waves 1-2 axe shown in Fig. 2. 
Three case examples were selected for study here, and are indicated by El, E2, and E3; 
each is 60 days in length and starting dates are 1 July 1981, 15 July 1984, and 1 March 
1987, respectively. Monthly mean 30 hPa QBO zonal winds reveal an apparent out of 
phase relationship with the eastward propagating components. A recent study by Shiotani 
and Horinouchi (1993) found a similar relationship, but between 70 hPa easterly zonal 
QBO winds and lower stratospheric Kelvin wave activity. The similarity suggests that 
eastward waves 1-2 in Fig. 2 are likely due to lower stratospheric Kelvin waves. This is 
supported further in the results by Randel and Gille (1991) which found a characteristic 
semiannual variability of Kelvin waves in the upper stratosphere and lower mesosphere with 
maximum activity near solstice. In contrast, semiannual variation of spectral amplitudes 
is not apparent in Fig. 2. 

Kelvin wave signatures in total ozone are expected to exhibit eastward propagating 
waves 1 or 2 with maximum amplitudes along the equator, and periods between 4 and 15 
days. The spectra in Fig. 3a show eastward waves 1-2 dominating all three episodes along 
the equator, with periods ~5-15 days. The line of constant phase speed shown in each 
frame is centered through the zonal wavenumber contributing most to total amplitude. 
Figure 3b shows structures typical of Kelvin waves, with eastward components exhibiting 
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strong meridional symmetry with a relative maximum along the equator. All three cases 
could be lower stratospheric “slow” Kelvin waves since phase speeds appear to be less than 
~40 ms 1 2 (see, for example, Wallace and Kousky, 1968). Faster phase speeds (exceeding 
60 ms -1 ) are more generally found in the upper stratosphere (Hirota, 1978; Salby et al., 
1984). 

Simple zonal wind advection of ozone does not explain these results. Monthly mean 
zonal winds at Singapore (not shown) for pressure levels 10, 15, 20, 40, 50, and 70 hPa 
indicate that westerly winds occur only at the lowest and highest altitudes for episodes 
El and E3, respectively, with values less than 12 ms -1 . Even at 12 ms -1 , the period of 
eastward propagating zonal wave 1 (2) would be 38 (19) days, which is roughly twice the 
observed values shown in Fig. 3. 

Figure 4 shows Hovmoller diagrams for episodes E1-E3. Bandpass filters were applied 
to retain the primary equatorial spectral responses shown in Fig. 3 (no eastward/ westward 
filtering is employed). Oscillation periods are ~6-10 days for each episode when amplitudes 
appear largest (exceeding 3 DU). Note eastward phase movement of strongest amplitudes, 
consistent with Kelvin wave dynamics. 

It is a possible that tropical clouds accompanying Kelvin waves could induce some 
of the features found here because roughly 10% of total column ozone is thought to lie 
below the tropopause and TOMS retrievals depend on the presence or absence of clouds. 
However, the model calculations presented above suggest that stratospheric ozone alone 
could account for the perturbation amplitudes observed. 

Eastward bandpass-filtered root-mean-square amplitudes (combined zonal waves 1-2) 
calculated for the three events are shown in Fig. 5. In each case an equatorial maximum 
is evident, with largest amplitudes for E2 and E3 near the African continent, and for El, 
along the western Pacific Ocean region. Longitudinal differences between episodes E1-E3 
could be caused by different tropical forcing regions, the effects of clouds and stratospheric 
winds, background gradients in ozone, or a combination of these. 


Summary 

1. Tropical Kelvin wave signatures in total column ozone are found in three separate 
case studies. Observed periods for zonal waves 1-2 were ~5-15 days, and amplitudes 
exceeded 3, and sometimes 5 DU. While tropical clouds associated with Kelvin waves 
could contribute to the observed amplitudes, a model calculation (independent of any 
cloud effects) shows that lower stratospheric ozone perturbations alone are sufficient to 
induce the observed TOMS Kelvin wave amplitudes. 

2. A strong relationship is found between the 30 hPa QBO wind and equatorial 
eastward zonal waves 1-2 in toted ozone, in contrast to semiannual variations in upper 
stratospheric Kelvin ozone waves found by earlier investigators. This implies that eastward 
Kelvin wave features in TOMS data sire influenced more by Kelvin waves in the lower-to- 
middle stratosphere than at higher altitudes. Furthermore, a previous study using 15 
years of tropical rawinsonde data showed that maximum Kelvin wave activity in the lower 
stratosphere occurs when the QBO zonal wind easterlies at 70 hPa are maximum. 
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FIGURE 1. Peak amplitudes in modeled zonal wave 1 (wave 2 results axe the same) 
Kelvin wave perturbation amplitudes (DU) in total column ozone plotted vs. vertical 
wavelength (bottom; units km) and intrinsic (u = 0) period (top; units days). UMKEHR 
levels 3-6, 7-10, and 3-10 refer to 5.3-127 hPa, 0-5.3 hPa, and 0-127 hPa, respectively. 
The latter is included for comparing cancellation effect(s) between the upper and lower 
stratosphere. 
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FIGURE 2. Calculated westward/eastward TOMS power spectra covering 1979-1991 
for combined zonal waves 1-2 along the equator, and 30 hPa monthly mean zonal winds 
(units ms -1 ) from Singapore (middle). A 180-day window was stepped forward in time 
at one-month intervals to compute spectra. Only the largest frequencies are shown, with 
periods in days indicated (15, 8, 5, and 3 days). Solid (dashed) contours in spectra are 
0.02,0.04,0.06,... (0.01) with 0.02 darkest; units are DU 2 -day. Symbols E1-E3 identify 
three 60-day case studies discussed in text. 
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amplitude of combined zonal waves 1-2. 
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FIGURE 4. Time vs. longitude Hovmoller diagrams for case studies E1-E3. Spatial 
filtering: zonal waves 1-2. El and E3 (E2) used a bandpass filter response with half- 
amplitudes at 15 and 4-day (20 and 4-day) periods. Solid (dashed) contours begin at zero 
(-1) and increment (decrement) by 1. Units: DU. 
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FIGURE 5. Latitude vs. longitude eastward root-mean-square amplitude plots for 
case studies E1-E3. Amplitudes are calculated from the eastward components of the data 
sets shown in Fig. 4. Solid (dashed) contours begin at 1.0 (0.9) and (increment) decrement 
by 0.1. Units: DU. 





